LAPTH-986 
KEK-CP-141 
KEK-TH-892 

Full 0(a) electroweak and 0{as) corrections to 

e+e- tiH 

G. Belanger^\ F. Boudjema^\ J. Fujimoto^-*, T. Ishikawa^-*, 
T. Kaneko^), K. Kato^), Y. Shimizu^), Y. Yasui^) 

1) LAPTH\ B.P.llO, Annecy-le-Vieux F-74941, France. 
2) KEK, Oho 1-1, Tsukuha, Ibaraki 305-0801, Japan. 
3) Kogakuin University, Nishi-Shinjuku 1-24, Shinjuku, Tokyo 163-8677, Japan. 



Abstract 

We present the full 0{a) electroweak radiative corrections to associated Higgs 
top pair production in e~^e~ collisions. We combine these results with a new cal- 
culation of the full one-loop QCD corrections. The computation is performed with 
the help of GRACE-loop. We find that the 0{a) correction can be larger than the 
0{as) corrections around the peak of the cross section especially for a light Higgs 
mass. At threshold these corrections are swamped by the QCD corrections which 
are enhanced by the gluon Coulomb contribution. We have also subtracted the 
complete QED corrections and expressed the genuine weak correction both in the 
a-scheme and the G^-scheme. This reveals that the genuine weak corrections are 
not negligible and should be taken into account for a precision measurement of this 
cross section and the extraction of the Yukawa ttH coupling. 
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1 Introduction 



After the discovery of the Higgs particle at the Large Hadron Colhder (LHC), one of 
the most pressing issue is a proper determination of the properties of this scalar since 
this would be an important window on the mechanism of electroweak symmetry breaking 
and the generation of mass. The LHC will be able to furnish a few measurements on the 
couplings of the Higgs to fermions and gauge bosonsfl but the most precise measurements 
will be performed in the clean environment of a future e~^e~ linear collider(LC) j2l El El- 
For example, from the measurement of the Higgs decay branching ratios, the Yukawa 
couplings of the light fermions can be determined at the per-cent level at a ^/s = 300 — 500 
GeV linear collider |2l El E] if the Higgs boson has a mass below the W pair threshold. 
This mass range for the Higgs is consistent with the latest indirect precision data [HI and 
covers the range predicted for the lightest Higgs of the minimal supersymmetric model 
(MSSM). At a TeV scale LC, the associate production of a Higgs boson with a top quark 
pair, e"*"e^ ttH, provides a direct information on the top-Higgs Yukawa coupling. In 
the Standard Model, SAi, the cross section of the e+e^ tiH process reaches a few 
fb, for a light Higgs and for centre of mass energies ranging from 700 GeV to ITeV. The 
expected accuracy for the determination of the top-Higgs coupling is of order 5% through 
the precision measurement of this process at the LC experiment [3 El El E] • Considering 
such a high accuracy on the ttH coupling one needs, on the theoretical side, to take into 
account the effect of radiative corrections. The purpose of this letter is to provide the full 
one-loop electroweak and QCD corrections to e^e~ — > tiH for a standard model Higgs. 
Preliminary results have already been presented in [7|. 

The full tree-level calculation of e^e~ tiH has been done a decade agojHj. An earlier 
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approximate calculation had been performed by only taking into account the (dominant) 
photon exchange diagrams[9j. QCD radiative corrections have also been performed by 
two groups. Dawson and Reina investigated the 0{as) corrections but only to the dom- 
inant photon exchange contribution JU]. The full 0{as) correction has been computed 
by Dittmaier et a/. Recently the supersymmetric QCD corrections have also been 
discussed in jT2]- On the other hand, due to the presence of the large top Yukawa cou- 
pling, the electroweak radiative corrections may also be sizable. However, the calculation 
of the electroweak radiative correction has been missing. We will present new results of 
the full 0{a) corrections consisting of virtual and soft corrections as well as hard photon 
radiation for the process e'^e'^ tiH in the SM. and will combine this result with the 
QCD corrections. 

2 Grace-Loop and the calculation of e^e~ tiH 

Our computation is performed with the help of GRACE-loop. This is a code for the 
automatic generation and calculation of the full one-loop electroweak radiative correc- 
tions in the SM.. It has been successfully tested for a variety of one-loop 2 — > 2 elec- 
troweak processes JHl- It also provided the first results on the full one-loop radiative 
corrections to e^e~ — » vvH ^1 which have recently been confirmed by an indepen- 
dent calculation [16J. For all electroweak processes we adopt the on-shell renormalisation 
scheme according to^Sj^JEl- For each process some stringent consistency checks are 
performed. The results, for the part pertaining to the electroweak corrections, are checked 
by performing three kinds of tests at some random points in phase space. For these tests 
to be passed one works in quadruple precision. Details of how these tests are performed 
are given in|18[ITlj. Here we only describe the main features of these tests. We first check 
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the ultraviolet finiteness of the results. This test applies to the whole set of the virtual 
one-loop diagrams. In order to conduct this test we regularise any infrared divergence by 
giving the photon a fictitious mass (we set this at A = lO^^^GeV). In the intermediate 
step of the symbolic calculation dealing with loop integrals (in n-dimension) , we extract 
the regulator constant Cuv = 1/^ — 7_e + log47r, n = A — 2e and treat this as a parame- 
ter. The ultraviolet finiteness test is performed by varying the dimensional regularisation 
parameter Cjjv- This parameter could then be set to in further computation. The 
test on the infrared finiteness is performed by including both loop and bremsstrahlung 
contributions and checking that there is no dependence on the fictitious photon mass A. 
An additional stability test concerns the bremsstrahlung part. It relates to the indepen- 
dence in the parameter which is a soft photon cut parameter that separates soft photon 
radiation and the hard photon performed by the Monte-Carlo integration. A crucial test 
concerns the gauge parameter independence of the results. Gauge parameter indepen- 
dence of the result is performed through a set of five gauge fixing parameters. For the 
latter a generalised non-linear gauge fixing conditional^] has been chosen. 



Cgf = -^\{d^ - ^eaA^ - igcwpZ^)W^+ + + SH + zkx3)x^\^ 



-id.Z + + ~eH)x,Y - ^{d.AY . (1) 

The X represent the Goldstone. We take the 't Hooft-Feynman gauge with = = 
= 1 so that no "longitudinal" term in the gauge propagators contributes. Not only 
this makes the expressions much simpler and avoids unnecessary large cancelations, but 
it also avoids the need for high tensor structures in the loop integrals. The use of the five 
parameters, a, /5, 5, k, e is not redundant as often these parameters check complementary 



sets of diagrams. Let us also point out that when performing this check we keep the full set 
of diagrams including couplings of the Goldstone and Higgs to the electron for example, 
as will be done for the process under consideration. Only at the stage of integrating over 
the phase space do we switch these negligible contributions. 

Although the system is not fully adapted for the computation of generic QCD correc- 
tions, it is quite straightforward to implement the QCD (final state) radiative corrections 
to e^e~ tiH . Indeed these corrections are rather QED-like corrections. The infrared 
divergence can be treated by giving the gluon an infinitesimal mass while the ultraviolet 
divergences are treated via dimensional regularisation. Also here we adopt an on-shell 
scheme in particular for the top mass and wave function renormalisation. The QCD tiH 
counterterm is then defined in terms of the top mass counterterm and the wave function 
constant. We have checked the infrared and ultraviolet finiteness of the QCD part also. 

The full set of the Feynman diagrams within the non-linear gauge fixing condition con- 
sists of 12 tree-level diagrams and 2327 one-loop diagrams (with 164 pentagon diagrams) 
for the electroweak 0{a) correction to the process e^e^ — > tiH, see Fig. ^for a selection 
of these diagrams. Even though we neglect the electron-Higgs coupling, the set of dia- 
grams still includes 6 tree-level diagrams and 758 one-loop diagrams (with 29 pentagon 
diagrams). We define this set as the production set. To obtain the results of the total 
cross sections, we use this production set. The handling of both the scalar and tensor 
pentagon integrals is done exactly along the lines in as was used for the calculation 
of e^e~ vvH . 

Our input parameters for the calculation of e^e~ tiH are the following. We will 
start by presenting the results of the electroweak corrections in terms of the fine structure 
constant in the Thomson limit with = 137.0359895 and the Z mass Mz = 91.187 GeV. 
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Figure 1: A small selection of different classes of loop diagrams contributing to e+e~ — >■ ttH . 
We keep the same graph numbering as that produced by the system. The first class of diagrams 
(the first 5 diagrams) are QED corrections, with the first row consisting of final state corrections. 
The QCD corrections diagrams can be inferred from the latter. The pentagon is a QED initial- 
final state interference. The second class groups genuine electroweak corrections including self- 
energy, triangle, box and pentagon corrections. Note that for the self-energy diagrams we collect 
all contributions including the counterterms in the blob. Note that the top self-energy Graph 
2260 also contains a photonic correction that should be included in the first class of diagrams. 
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The on-shell renormalization program uses My/ as an input. However, the numerical value 
of Mw is derived through Ar[in] with = 1.16639 x lO^^GeV"^*. Thus, Mw changes as 
a function of M^- For the lepton masses we take rrie = 0.510999 MeV, = 105.658389 
MeV and rrir = 1.7771 GeV. For the quark masses beside the top mass Mt = 174 GeV, we 
take the set = Ma = 63 MeV, M, = 92 MeV, = 1.5 GeV and Mb = 4.7 GeV. With 
this we find, for example, that Mw = 80.3759 GeV for Mh = 120 GeV and Mw = 80.3469 
GeV for Mh = 180 GeV. For the QCD coupling, we choose as{Mz) = 0.118 as an input 
and evaluate as{Mt) = 0.10754 with the next-to-next-to-leading order renormalization 
group equation. 

As well known, from the direct experimental search of the Higgs boson at LEP2, 
the lower bound of the SAi Higgs boson mass is 114.4 GeV[2I]. On the other hand, 
indirect study of the electroweak precision measurement suggests that the upper bound 
of the SAi Higgs mass is about 200 GeVj3]. In this paper, we therefore only consider a 
relatively light SAi Higgs boson and take the two illustrative values Mh = 120GeV and 
Mh = 180GeV. 

Let us first present some quantitative consistency tests on our results. For the elec- 
troweak part, the ultraviolet finiteness test gives a result that is stable over 20 digits when 
one varies the dimensional regularisation parameter Cuv As for the gauge parameter in- 
dependence checks, our results are stable over 26 digits when varying any of the non-linear 
gauge fixing parameters. For the QED infrared finiteness test we also find results that 
are stable over 20 digits when varying the fictitious photon mass A. As for the kc stability 

test our results are consistent within a Monte-Carlo statistical error of 0.02%. 

*The routine we use to calculate Ar has been slightly modified from the one used in our previous 
paper on e~^e~ — *■ vDH |14j to take into account the new theoretical improvements. It reproduces quite 
nicely the approximate formula in |20|. 
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We also checked the tests of the stabihty on the ultraviolet and the infrared finiteness 
for the QCD calculation. The ultraviolet finiteness test gives a result that is stable over 
30 digits. The sum of loop and bremsstrahlung contributions is stable over 13 digits when 
varying A^, the infrared gluon mass regulator, independence is consistent with a Monte- 
Carlo statistical error of 0.2%. In addition, we reproduced the previous results of ref. [TUj 
when setting to zero the Z exchange diagrams and also exactly reproduced the results 
of Dittmaier et a/.^T] with the full 0{as) calculation within the on-shell renormalization 
scheme. 

3 Results 

3.1 Full 0(a) and 0{as) corrections 

At tree-level, for the Higgs masses we are considering, the cross section shows a steep 
rise just after threshold and slowly decreases past the maximum of the cross section, 
see Fig. Efa). For the measurement of the tiH coupling it is most useful to run at the 
maximum of the cross section. For example for = 120GeV this maximum occurs 
around ^/s = 700 — SOOGeV where the cross section is in excess of 2/6, see Fig. El^a). For 
a total integrated luminosity of lab~^ the Icr statistical error corresponds to about a 2% 
precision. Thus the theoretical knowledge of the cross section at 0.2% is quite sufficient. 
It is important to keep in mind that the dominant contribution to e~^e^ tiH for 
the energies we are considering is due to the photon exchange diagram, the Z exchange 
diagrams with Higgs radiation off the Z or the top is much smaller. This is especially true 
at threshold. For example for Mh = 120GeV and ^/s = 500GeV, the photon exchange 
diagram alone contributes 90% of the total cross section whereas Higgs radiation off the Z 
is less than 0.2%. The importance of the photon exchange contribution lessens somehow 
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Figure 2: (a) Total cross section as a function of the centre of mass energy for Mh = 
120 GeK We show the total cross sections for the tree level, full 0{a) and 0{as) level in 
(a). The relative corrections are shown in (h). Solid lines are tree level, dashed lines are 
the full 0{a) and dotted lines are the 0{as) corrections. We take as{Mt) = 0.10754. In 
addition, the genuine weak correction 5w o,nd the relative correction 5y^ in the scheme 
are presented. 
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at high energies, but even at ITeV this contribution accounts for about 80% of the total. 

The full 0{a) electroweak correction and the full 0{as) QCD inclusive cross section 
take into account the full one-loop virtual corrections as well as the soft and hard bremm- 
strahlung contributions. These corrections are shown Fig. |2^a). The relative corrections 
defined as 

<JOia,a^) — -L 

^tree 

are shown in Figure |21^b). 

Let us first turn to the inclusive QCD correction. It is useful to write this correction 

as 

So{a.) = Cf A„ Cf = 4/3. (2) 

The QCD correction has a scale dependence which however, at this order, is fully contained 
in a<j(/i). We will provide the exact value of ^^(/i) together with So(as) so that a comparison 
with other calculations is straightforward. As already noted in [TTH ITT] we confirm that 
the QCD corrections are quite large at threshold increasing the tree-level cross section by 
about 50%, for Mh = 120GeV and y/s = SOOGeV. This large increase is due essentially to 
the gluonic Coulomb corrections. At energies where the cross section is at its highest and 
the process is most likely to be of interest, the QCD corrections are quite modest. For 
example, for Mh = 120GeV and ^/s = SOOGeV, as seen also in Tabled the correction is 
a mere ~ —1%, the residual scale dependence is very small. Indeed as shown in Tabled 
the corrections for a scale /i = Mf and /i = i/i = 800GeV~ 4Mt, vary from —1. to 
—0.8%. A similar trend is observed for Mh = 180GeV, see Table ^ As the energy 
increases further past the maximum of the cross section, these corrections turn negative 
but remain moderate at around —5% for ^/s = ITeV. At these energies for the light Higgs 
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masses we are considering the Mh dependence of the corrections lessens. 





Mj^(GeV) 


entree (fb) 


o"ci(Q4A/t))(fb) 5o{ 


as{Mt)){%) 




600 GeV 


120 


1.7293 ±0.0003 


1.977 ±0.001 


14.3 


12.4 




180 


0.33714 ±0.00004 


0.4383 ±0.0004 


30.0 


26.0 


800 GeV 


120 


2.2724 ±0.0005 


2.250 ±0.002 


-1.0 


-0.8 




180 


1.0672 ±0.0003 


1.0856 ±0.0007 


1.7 


1.4 


1 TeV 


120 


1.9273 ±0.0005 


1.812 ±0.003 


-6.0 


-4.9 




180 


1.1040 ±0.0003 


1.049 ±0.001 


-5.0 


-4.1 



Table 1: QCD corrections for e^e ttH . We also display the Monte-Carlo integration 
errors. We consider two different schemes. (A) We choose the renormalization scale /i 
of the QCD coupling a, at Mt = 174 CeV with as{Mt) = 0.10754. (B) We take = y/s 
with as = 0.09330 at ^ = 600 CeV, as = 0.09051 at ^ = 800 CeV and as = 0.08847 
at ^ = ITeV 



We now turn to the total 0{a) electroweak corrections. Again taking as an example 
Mh = 120GeV, the 0{a) electroweak correction is about —7.5% around threshold at 
y/s = 500 GeV and is therefore swamped by the QCD correction. However as the energy 
increases, contrary to the QCD correction, the full electroweak correction slowly increases 
and turns positive for ^/s > 600 GeV. Around the maximum of the cross section at 
^/s = 800, the full electroweak correction is about ±4% and thus larger than the full 
QCD correction (—1%). For yet higher energies these corrections tend to cancel each 
other. For example for M^ = 120GeV this occurs around y/s = ITeV. We also show 
some numerical results Mh = 180 GeV in Table El The 0{a) corrections for the Higgs 
mass of 180GeV are rather small in the range ^/s = 800 GeV — 1 TeV, slightly increasing 

from 2% to -0.5%. In particular, for = ITeV and Mh = 180GeV the full 0{a), 

because of its smallness, may not be resolved from the corresponding QCD corrections 
due to the scale dependence of the latter. 
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C^iree (fb) 


CrO(a) (fb) 




600 GeV 


120 


1.7293 ± 0.0003 


1.738 ±0.002 


0.5 




180 


0.33714 ±0.00004 


0.3126 ±0.0003 


-7.3 


800 GeV 


120 


2.2724 ±0.0005 


2.362 ±0.004 


3.9 




180 


1.0672 ±0.0003 


1.050 ±0.002 


-1.6 


1 TeV 


120 


1.9273 ± 0.0005 


2.027 ±0.004 


5.2 




180 


1.1040 ±0.0003 


1.098 ±0.002 


-0.5 



Table 2: As in the previous table but for the total 0{a) electroweak corrections. 
3.2 The genuine weak correction 

In order to quantify the genuine weak corrections one needs to subtract the full QED 
corrections from the full 0{a) corrections. This is important because it is well known 
that the QED corrections can be quite large and that in e+e^ processes those from the 
initial state need to be resummed|22j. For the process at hand, which at tree-level proceeds 
through s-channel neutral vector bosons, these QED corrections form a gauge invariant 
set. This set may be further subdivided into three subsets that are also separately gauge 
invariant: i) initial state radiation, ii) purely final state radiation and Hi) the initial-final 
state QED interference. 

i) The dominant initial state QED virtual and soft bremsstrahlung corrections are given 
by the universal soft photon factor that leads to a relative correction 14] 

= ^ [{Le - 1) In |- ± ^L. ± ^ - l) , = ln(./m^) . (3) 

where me is the electron mass, the beam energy (s = 4:E'^) and kc is the cut on the 
soft photon energy. 

ii) The total QED final state radiation can also be read off from the result of the QCD 
radiative correction through the replacement as{fi) Cp ^ a ui Eq. |21((5t is the electric 
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charge of the top). 

in) The initial-final state QED correction is ultraviolet finite. Within our system this 
contribution can be easily isolated and combined with the appropriate bremmstrahlung 
counterpart leading to an infrared finite result. 

Although this approach of extracting the full QED correction is the most simple one, 
we have also calculated the full QED corrections separately and subtracted their con- 
tributions from the full 0{a). In order to perform this subtraction, the QED virtual 
corrections are generated by dressing the tree-level diagrams with one-loop photons (the 
photon self-energy is not included in this class). Moreover one needs to include some 
counterterms. One only has to take into account the purely photonic contribution to the 
top mass counterterm as well as the wave function renormalisation constants of the elec- 
tron and the top. Performing this more direct computation, we confirmed, that especially 
around threshold, to a large extent the bulk of the QED corrections originate from the 
initial state universal corrections. Moreover this also checked the extraction of the final 
QED corrections. A break up of the soft and virtual QED corrections into initial, final 
and interference is shown in Table El 
We define the genuine weak relative correction as, 

r X xQED r s:QED irQED 

OW = do{a) - d = do{a) " (Jy+S - (Jftarrf • 

^v+s is complete QED virtual and soft correction whereas S'j^^jj' is the hard photon 
contribution. The weak corrections after subtraction of the QED corrections are shown 
in Fig. 121 for = 120GeV. The total genuine weak corrections are not small, being 
largest around threshold (~ +25% at y/s = SOOGeV) and decrease monotonically as the 
energy increases, see Fig. |2fb). At = ITeV they reach about 6%. These corrections 
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Mh 
(GcV) 


(GcV) 


i'^hard^Fulli'^)) 


QEU 
^V+S.Init. 

(fb) 


QEU 
^V+S.Fin. 

(fb) 


QEU 
^V+S.Int. 

(fb) 


§QED 
(%) 


(%) 


120 


















600 


-2.6092 
(2.3333) 


-2.557 


-0.012 


-0.043 


-16.0 


16.5 




800 


-3.6667 
(3.5391) 


-3.516 


-0.055 


-0.099 


-5.6 


9.5 




1000 


-3.2622 
( 3.2507) 


-3.086 


-0.071 


-0.109 


-0.6 


5.8 


180 


















600 


-0.50490 
( 0.41839) 


-0.4985 


0.0056 


-0.0072 


-25.7 


18.4 




800 


-1.7120 

( 1.5975) 


-1.651 


-0.020 


-0.043 


-10.7 


9.1 




1000 


-1.8589 
(1.8048) 


-1.768 


-0.035 


-0.058 


-4.9 


4.4 



Table 3: Extraction of the QED corrections. CFy^g p^n corresponds to the cross section 
for the full one-loop QED virtual and soft bremsstrahlung with kc — 0.001 GeF. crv+s,init. 
extracts the initial state radiation, cry^g p^^^ gives the final state QED correction whereas 

Cv+sint initial-final QED interference contribution. We also give cr'^ardFuU which 

is the full hard photon radiation cross section. All cross sections are in fb. We also 
give the relative QED correction (after including hard radiation) as well as the relative 
genuine weak correction as defined in the text. Note also that the extraction of the total 
^v+s,Fuii ^"i^^ ^hard,Fuii been computed with higher accuracy than the individual ( S-h V) 
contributions. We can check that the two computations (full) and adding the individual 
contributions agree. 



13 



could therefore, for this Higgs mass, always be disentangled from the QCD corrections. 
Past y/s = 600GeV and up to y/s = ITeV , they are larger, in absolute terms, than 
the QCD corrections. A similar trend also occurs for Mh = ISOGeV, see Table El For 
^/s = 600GeV we find 6w ~ +18% this correction drops with energy reaching Sw ~ +4% 
at ITeV where it almost cancels the corresponding QCD correction. 

Having subtracted the genuine weak corrections one could also express the corrections 
in the scheme by further extracting the rather large universal weak corrections that 
affect two-point functions through Ar. This defines the genuine weak corrections in the 
G^j scheme as 6yy = 5iy — 3Ar. For reference, one has Ar = 2.55% for = 120 GeV and 
Ar = 2.70% for Mh = 180 GeV. For e^e~ — > uuH this procedure helps absorb a large 
part of the weak corrections. Another advantage is that much of the (large) dependence 
due to the light fermions masses also drops out. For Mh = 120GeV, the relative correction 
6^ is shown as a function of energy in Fig. |21 Adopting this scheme, we find that in 
the energy range where the cross section is largest, ^/s = TOOGeV to ^/s = ITeV, the 
correction remains modest changing from about 5% to —2%. These corrections could 
be "measured above" the QCD corrections. For Mh = 180GeV, the correction ranges 
from 1% to —4% as the energy changes from 800GeV to ITeV. However at energies 
around threshold, the genuine weak corrections in the scheme are large (and positive), 
although about a factor 3 to 2 smaller than the QCD corrections. For Mh = 120 GeV the 
genuine non-QED correction in the scheme reaches +17% at y/s = 500GeV (~ SOGeV 
above threshold). They are about 10% for Mh = 180GeV = 600GeV (~ 70GeV above 
threshold). These corrections slightly decrease with energy following a trend similar to the 
QCD corrections although the decrease is not as fast. One knows from e'^e'^ tt [^ 1^ 
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that there is a Yukawa counterpart to the Coulomb large gluonic correction mediated 
by a Higgs exchange which is large at threshold for a large top mass and a very light 
Higgs mass. This phenomenon can only partially account for the large 17% increase of 
the cross section at ^/s = SOOGeV for Mh = 120GeV, considering the value of the Higgs 
mass. We would like to argue that because the tree-level cross section is dominated by 
photon exchange, a description of the cross section in terms of instead of a may 
not be the most appropriate. Rather, the photon couplings should be parameterized in 
terms of the running a. Therefore especially at threshold one should subtract from 6w 
2Aa{s)+Ar instead of 3Ar. For illustration we will only use the light fermion contribution 
to Aa{s). The correction defined this way, Smixed = ~ 2Aa;(s) — Ar, brings down the 
remaining weak contribution to Smixed ~ 7% for Mh = 120GeV at -^i = SOOGeV and 
1.4% for Mh = ISOGeV at = 600GeV. For higher energies, above = 700GeV, 
this prescription gives large (negative) corrections. For instance, for y/s = ITeV we find 

Smixed 12% for Mh = UOGeV and Smixed 14% for Mh = ISOGeV. Therefore none 

of the prescriptions, or Smixed reproduces the full genuine weak corrections across the 
whole energy range from threshold to ITeV. This should not be surprising as it has been 
known, already for 2-^2 processes, that the contribution of boxes becomes important as 
the energy increases and that the scheme is not always the most appropriate. A case 
in point is e~^e~ — > ZH [25j. One may also enquire whether the leading corrections to 
the tiH coupling could account for most of the weak radiative corrections. The leading 
corrections to the ffH vertex, in the G^ scheme, had been worked out |2nil2Zl- The 
corrected tiH vertex yuH writes in terms of the tree- level one, y^^n^ 




(4) 
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In the total cross section with Mt — 174GeV this only accounts for about 2.2% weak cor- 
rection. Therefore one sees that to properly take into account the electroweak corrections 
to the e+e" ttH a full calculation is needed. 

4 Conclusions 

We have performed a full one-loop correction to the process e+e" ttH which, at a 
future linear collider running in the energy range 700GeV to ITeV, can allow a direct 
determination of the important Yukawa coupling ttH. The full one-loop corrections com- 
bine both the full electroweak corrections as well as the QCD corrections. Apart from 
the ultra-violet and infrared finiteness tests of the results, we performed, for the more 
involved electroweak sector, an extensive gauge parameter independence check, kc sta- 
bility has also been verified, kc is the photon energy parameter that separates soft and 
hard photon (and gluon) radiation. For the electroweak part we have also extracted the 
contribution of the QED corrections. The final state QED correction has also been used 
as a further check on the QCD part. This extraction helps define the genuine weak cor- 
rections. We have expressed the latter both in the a-scheme and the G^-scheme. We 
find that, for all the Higgs masses that we have studied, the full 0{a) corrections are 
swamped by the large QCD corrections at threshold. The latter can increase the cross 
section by as much as 50% due to the threshold Coulomb enhancement. It is interesting 
to note that after the extraction of the QED corrections, the genuine weak corrections 
at threshold are important (~ 20%), although still small compared to the corresponding 
QCD corrections. However, around the peak of the cross section where this process is 
most likely to be of interest for the extraction of the Yukawa coupling, the electroweak 
corrections can dominate over the QCD corrections. In this energy range both corrections 

16 



are under control and for Mh = 120GeV, the genuine weak corrections expressed in the 
Gfj, scheme are modest ranging between 5% to —2% for ^/s = SOOGeV to = ITeV. 
Nonetheless they have to be taken into account for a precision analysis of this process at 
a future linear collider. 
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Addendum 

Our results for the full 0{a) corrections were first reported in C|. In the meantime while 
finalising this article there appeared a paper on the same subject ||28j. These authors 
calculate the 0{a) corrections using the system FeynArts and FeynCalc|29j but do not 
attempt to isolate the genuine weak corrections. These authors use a rather different set 
of input parameters. For the light quark masses they use the current quark masses which 
would give a much too large value of Aa(M|). Nonetheless taking the same input pa- 
rameters we have run our program for a few points reported in|2H]- With = llSGeV, 
for both ^/s = SOOGeV and ^/s = lOOOGeV we find very satisfactory agreement, at least 
as far as one can read from their graph (Fig. 5 of [^H"). We obtain — 4.2%(y^ = SOOGeV) 
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and 7.3% (-yi = ITeV). However for higher energies we have not been able to reproduce 
their results. For instance for i/i = 2TeV we find So(^a) ~ 8% whereas they find ~ 4%. 

Note added 

After submitting this paper there appeared another calculation of this processj^U]. This 
new calculation perfectly confirms all of our results, for instance the relative corrections 
in the two codes agree within 0.1% at all energies and for all masses. 
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